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基于 调制 宽带 转换 器 的 频谱 分 辨 率 研究 
EC MEA! Hy Aw? 
(1. 云南 大 学 信息 学 院 ， 云 南 昆明 650500; 2. 中 国 科 学 院 云南 天 文 台 射 电 天 文 与 VLBI1 研究 团 组 ， 
云南 昆明 650216) 
摘要 : 目前 ， 压 缩 感知 的 研究 热点 主要 集中 在 重 构 算 法 改进 方面 ， 即 如 何 降 低 重 构 误差 和 重 构 
时 间 ， 且 均 在 实际 应 用 中 所 允许 的 范围 内 ,而 对 于 信号 频 域 分 析 中 的 频谱 分 状 率 提升 技术 的 关注 较 


少 。 本 文采 用 调制 宽带 转换 器 对 稀 足 信号 进行 采样 ,采用 正 交 匹配 追踪 算法 进行 重 构 输 入 信号 ， 通 


过 延长 时 域 采样 时 长 ， 实 现 单个 通道 低速 采样 的 同时 提升 频谱 分 状 率 。 对 低频 射电 天 文 信号 进行 仿 
真实 验 结果 表明 ， 在 单个 通道 采样 速率 相同 时 ， 时 域 采 样 样本 数 由 459 增加 至 999 时 ， 频 谱 分 辨 率 
由 0. 4357MHz 提升 至 0. 2002MHz， 但 采样 系统 的 总 采样 样本 数 仅 由 102 增加 至 222。 本 文 的 研究 为 
射电 天 文 信 号 等 高 频 窜 带 信号 处 理 提 供 了 很 好 的 借鉴 。 

关键 字 : 调制 帘 带 转换 器 ， 频 谱 分 辨 率 ， 压 缩 感 知 ， 射 电 天 文 信号 


中 图 分 类 号 : TN911. 7 


0 引言 

频谱 分 析 是 信号 处 理 的 重要 技术 手段 之 一 , 通过 对 频谱 构成 和 分 布 进 行 分 析 ， 可 获取 信号 的 某 
些 关 键 特征 。 针 对 高 频 窄 带 信号 ， 需 在 满足 香农 采样 定理 的 要 求 下 ， 采 用 较 高 速率 ADC 对 信和 号 采样 
和 进行 样本 数据 分 析 处 理 。 目 前 ， 最 高 速率 的 ADC 无 法 满足 学 者 对 频谱 分 析 的 需求 ， 因 此 ， 国 内 外 
学 者 对 信号 的 结构 和 特点 进行 分 析 , 针对 自身 或 通过 数学 变换 后 具有 稀 下 性 质 的 信号 ， 提 出 了 压缩 
感知 (Compressed Sensing, CS) 理论 “”。 压 缩 感知 理论 主要 包含 前 端的 压缩 采样 系统 和 后 端的 
音 号 重 构 两 个 部 分 ， 相 比 传统 采样 技术 ， 压 缩 感 知 理论 将 采样 和 压缩 合 二 为 一 ， 降 低 了 采样 速率 、 
数据 传输 带宽 以 及 存储 空间 。 近 年 , 国内 外 的 学 者 对 压缩 感知 的 研究 内 容 为 重 构 算 法 的 改进 和 实际 
点 用 , 在 实际 应 用 的 误差 允许 范围 内 , 压缩 感知 相关 理论 能 以 较 低 的 采样 速率 对 稀疏 信号 作 亚 奈 奎 
斯 特 采样 ， 并 能 完美 的 重 构 出 输入 信号 。 但 针对 高 频 窄带 信号 处 理 的 压缩 感知 理论 指导 实际 应 用 的 
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研究 不 足 , 在 理论 研究 中 的 输入 信和 号 模型 是 时 域 长 度 设 定好 的 信号 序列 ， 既 频谱 的 分 辨 率 也 因此 设 
定好 了 ， 未 能 有 效 的 指导 实际 应 用 中 如 何 提升 信号 频谱 的 分 辨 率 。 

在 天 文 观测 中 , 射电 天 文 信号 的 谱 线 观测 成 为 研究 宇宙 的 重要 方式 , 通过 对 谱 线 进行 分 析 可 以 
诊断 有 关 天 体 的 基本 物理 条 件 ， 如 动能 温度 、 总 粒子 数 密度 、 速 度 、 磁 感应 强度 等 ”。 射 电 天 文 信 
号 的 谱 线 观测 对 终端 设备 的 频率 分 辩 率 也 提出 了 很 高 的 要 求 *， 为 了 提高 观测 谱 线 的 频谱 分 辨 率 ， 
需 在 采样 速率 一 定 的 情况 下 增加 时 域 的 采样 时 长 , 需 消 耗 处 理 器 内 部 的 乘法 器 、 累 加 器 以 及 存储 单 
元 等 资源 。 在 实际 的 射电 天 文 观 测 中 ， 仅 对 一 些 指定 的 谱 线 感 兴 些 指定 的 谱 线 在 频 域 上 具有 
稀疏 特性 ， 因 此 ， 本 文 把 压缩 感知 理论 应 用 到 射电 天 文 信号 观测 中 ， 通 过 对 前 端的 调制 宽带 转换 器 
(Modulated Wideband Converter, MWC)“ 采 样 系统 和 的 正 交 匹配 追踪 (0rthogonal Match Pursuit, 
OMPO " 重 构 算 法 进行 分 析 ， 实 现 低速 ADC 对 高 频 窄带 信号 的 压缩 采样 ， 并 可 接近 无 失真 地 重 构 出 
输入 信号 ， 仿 真实 验 表 明 ， 在 采样 点 数 低 于 奈 奎 斯 特 速 率 采样 样本 数 的 情况 下 ， 本 文 方法 有 效 地 提 
升 了 频谱 分 辨 紊 。 


1 MWC 采样 系统 和 重 构 算法 
1. 1 MWC 采样 系统 
MWC 是 由 Mishali 和 Eldar 提出 的 亚 奈 奎 斯 特 采样 系统 ， 并 对 系统 进行 了 详细 的 数学 推导 过 
程 , 其 由 多 个 通道 组 成 , 每 个 通道 由 伪 随 机 序列 发 生 器 、 低 通 滤波 器 和 低速 ADC 组 成 , 如 图 1 所 示 。 
MWC 采样 系统 主要 参数 有 : m 为 采样 系统 的 通道 数 ，T, 为 伪 随 机 序列 p(t) 的 周期 ，T 为 ADC 采 
样 周 期 ，aijx 为 第 i 通道 上 pi(t) 的 第 k 序 列 的 值 Cay € {41,-1}) ，M 为 pi(t) 一 个 周期 内 包含 +1 的 数 


由 


Hi, 低 通 滤波 器 h(t) = T7 sinc (Fe) Gill AREAS = 7 — 770. favo REM RAR. 


中 ，pi(t) 的 表达 式 为 : 


N 


pi(t) = ag, kT,/M € t < (k + 1)T,/M,0 <k <M —1 a) 


其 傅 里 叶 级 数 展开 表达 式 为 : 


—j221t 
pi(t) = pu oP 


j= ik Tp fol 


Cü — zu zu pi(t)e egt dt = RAE d aige Trp! “fme Trp! "dt (2) 


图 1 MWC 采样 系统 


Fig.1 The sampling system of MWC 


由 图 1 可 得 如 下 时 域 数学 表达 式 : 


Xi(t) = x(t)pi(t). yi(t)= (t) 8 xi(t). yiln] = yi(n7;) (3) 
TH, QUARERuSSE. SARIRO ey f BH. YP) Seas SB DBDE S8 HA E GR y; [n] I 


离散 时 间 傅 里 时 变换 CDTFTO ， 则 公式 (3) 的 频谱 表达 式 为 : 


Y()- X .xInle 27" —y | cuX(f — Up). f e[-5/2.5/2] (4) 


xut, Ly = ELEL] - 1, MERENTAS. Vm Att CA REA CU 
p 


y(f) = Az(f) (5) 


yf) = A YOA, Ym PI" 
z(f) = AG) Zh) 2 
zi(f) = X(f + (i — Lo + Dfy) 


A= [Cilmxn L =21)+1 
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1. 2 重 构 算 法 

MWC 采样 系统 的 后 端 重 构 过 程 是 先 计算 采 样 系统 各 通道 之 间 输 出 压缩 数据 的 协 方差 ,再 采用 仿 
禁 算法 找 出 包含 有 用 信息 的 频谱 片段 和 对 应 字典 的 索引 。 因 MWC 每 个 通道 上 输出 数据 的 长 度 是 有 限 
的 ， 令 数据 长 度 为 N， 则 各 通道 输出 数据 之 间 的 协 方差 矩阵 方程 为 : 


= 


Rie = 2n ZEN yilndys n] (6) 


Eldar 等 利用 特征 分 解 方法 把 对 z( 亡 频谱 片段 求解 过 程 转化 为 压缩 感知 问题 ， 即 对 有 进行 特征 


R= U*AU s (U*AV?)(u*AV?y = VV" 


V —AS (7) 


SUE, AZERIRHIEREZELSDSE AEE, UER RE EEE, OR APE TE PES 8 F AMEE 


HE OMO 问题 。 重 构 过 程 的 核心 步骤 为 通过 贪 梦 算 法 找 出 4 中 尽 可 能 少 的 列 向 量 和 $ 中 包含 有 用 信 
息 的 行 向 量 来 表达 出 Y， 把 少量 列 向 量 矩 阵 4s 的 伪 首 矩阵 (484s) -4148 带 入 公式 〈5) 中 进行 信号 重 
构 。 


2 频谱 分 辩 率 

随 着 芯片 制造 工艺 的 提升 ，ADC 采样 速率 可 达 几 个 甚至 十 几 个 GHZz， 信 号 处 理 的 总 带宽 得 到 很 
大 的 扩展 ， 但 对 高 频 罕 带 信号 的 谱 线 轮廓 结构 的 观测 依然 难以 满足 ， 比 如 ， 目 前 ADC 的 最 高 采样 速 
率 很 难 满 足 射电 天 文 观 测 中 高 频 窗 带 信 号 总 带宽 的 要 求 。 根 据 传 里 叶 变 换 原理 ,在 传统 奈 奎 斯 特 采 
样 方法 中 ， 采 样 时 间 分 辨 率 4t 与 接收 机 总 带宽 B 的 关系 为 4t. B = 1， 在 观测 时 间 r 内 的 采样 点 数 为 
N= 1/At =T:B。 根 据 高 斯 统计 ，N 个 采样 样本 的 总 误差 4T 是 单个 采样 的 误差 (系统 误差 ) Tsys 的 
1/VN， 接 收 机 的 极限 灵敏 度 为 4T/Tsys = 1/VzB “”， 提 高 接收 机 总 带宽 可 以 提升 极限 灵敏 度 ， 但 
意味 着 更 高 的 采样 速率 。 在 实际 射电 天 文 观测 中 ，ADC 的 采样 速率 不 小 于 奈奈 斯 特 速 率 ， 并 以 延长 
接收 机 接收 信号 的 时 长 来 提升 频谱 分 辨 率 , 增加 了 采样 样本 数 和 后 端 处 理 资源 的 消耗 。 较 高 的 频谱 
分 辨 率 ， 精 细 的 信号 频谱 结构 ， 对 分 析 和 研究 信号 特征 具有 重要 意义 。 


频谱 分 辨 率 4F、 采 样 频率 及 和 采样 样本 数 N 三 者 的 关系 为 : 


F, 


Af == (8) 


由 公式 CB) 可 知 ， 可 通过 降低 采样 速率 或 增加 采样 时 长 来 提高 频谱 分 辩 率 。 在 MWC 采样 系统 
中 ， 令 伪 随 机 序列 pi(t) 的 周期 为 M -Tuyo (Twyo 为 奈 夺 斯 特 采 样 周 期 , Tyyo = 1/jvro) ， 每 个 通道 


ADC 采样 速率 为 fyyo/M，MWC 的 总 采样 速率 降低 m/M , (m < M) 倍 ， 频 谱 分 辩 率 提升 M/m 倍 。 


3 实验 仿真 与 分 析 
3.1 基于 MWC 的 理论 分 析 


实验 仿真 信号 为 稀疏 多 频带 信号 ， 前 端 采 用 MWC 采样 系统 进行 采样 ,输出 压缩 数据 ， 后 端 采 用 
OMP 算法 对 压缩 数据 进行 运算 ， 重 构 出 输入 信号 。 稀 玻 多 频带 信号 的 参数 为 : 子 频带 数 N 为 8， 带 


宽 均 B 为 5MHz， 总 占 频 带宽 为 106Hz《〈 即 : 奈 奎 斯 特 采 样 频率 和 Avro) ， 信 品 比 为 20。MWC 采样 系统 


参数 为 : 通道 数 m 为 20 (m > 2N)， 伪 随机 序列 周期 长 度 M 为 195, f, = fp = ro/M ~ 51.28MHz， 


L = M = 195(Lo = 97). Y MWC 单个 通道 上 采样 样本 数 为 K ,公式 (4) 中 y(f) 和 矩阵 大 小 为 m x K, A58 
阵 大 小 为 m x M，z(f) 和 矩阵 大 小 为 M XK, yf) = 4z(f)。 输 入 信号 时 域 采 样 样本 数 和 频谱 分 辨 率 
关系 如 表 1 所 示 。 


表 1 时 域 采样 样本 数 和 频谱 分 辨 率 关 系 


Tab.1 The relationship between time domain sample number and spectral resolution 


单个 通道 采 
. 101 201 301 401 501 601 701 801 901 1001 
样 样本 数 
等 效 总 采样 
19695 39195 58695 78195 97695 117195 136695 156195 175695 195195 
样本 数 
频谱 分 辨 率 
507.7 255.1 170. 4 127.9 102. 4 85.3 73: 2 64.0 56.9 51.2 
/KHz 
重 构 误差 0.0029 0. 0023 0. 0014 0. 002 0. 0008 0. 0005 0. 0005 0. 0004 0. 0003 0. 0003 


MWC 采 样 系统 中 伪 随 机 序列 的 频率 要 不 小 于 奈 达 斯 特 采样 速率 才能 近乎 无 失真 的 重 构 出 输入 信 
号 ， 伪 随机 序列 在 硬件 实现 技术 上 的 难度 远 低 于 ADC， 因 此 在 采样 系统 中 采用 高 速 伪 随 机 序列 生成 
器 对 输入 信号 做 数字 离散 化 ， 经 低 通 滤波 器 后 进行 低速 采样 。 在 仿真 中 ， 伪 随机 序列 的 频率 等 于 对 
奈 斯 特 采 样 速率 ，MWC 采样 系统 总 采样 样本 数 为 mK， 压 缩 数 据 经 数学 变换 后 的 z( 放 矩阵 样本 数 为 


MK, Jl] MWC 采样 系统 的 等 效 总 采样 样本 数 为 MK。 因 此 ， 从 表 1 得 出 MWC 采样 系统 的 等 效 总 采样 


道 采 样 样本 数 的 增加 ， 


E 


样本 数 是 单个 通道 采样 样本 数 的 M 倍 , BITRE EAM as 随 着 单个 
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等 效 总 采样 样本 数 和 频谱 分 辨 率 随 之 增加 。 


3. 2 基于 MWC 的 低频 射电 天 文 信号 仿真 

低频 射电 天 文 信号 的 数据 来 源 于 云南 天 文 台 ， 以 200MHz 采样 速率 对 频段 在 55Mhz-65Mhz 的 射 
天 文 信号 进行 采集 的 多 组 数据 ， 本 实验 随机 选取 其 中 的 一 组 时 域 采样 数据 进行 仿真 ， 其 采样 点 数 
73 1024, 频谱 如 图 2 中 所 示 。MWC 采样 系统 的 参数 设置 为 : 通道 数 m 为 6, 伪 随 机 序列 周期 长 度 M 为 
27, f, = fo = fuyg/ M = 7.4074MHZ, L=M 227 (Lo 2 13). ， 总 采样 速率 为 44. 4444MHZ。 本 次 
实验 中 的 MWC 非 全 言 采样 系统 ， 需 要 设置 输入 信号 的 子 频带 数 为 2， 并 设 定 MWC 单个 通道 上 采样 样 
PK 为 不 同 值 时 的 仿真 数据 及 频谱 图 ， 如 表 2 和 图 2 所 示 。 

表 2 是 K 分 别 为 17、27、37 时 ， 对 应 MWC 采样 系统 采样 样本 数 、 等 效 总 采样 样本 数 、 频 谱 分 
辩 率 、 重 构 误 差 和 OMP 算法 的 重 构 时 间 的 仿真 实验 数据 。 通 过 表 2 可 以 看 出 ， 时 域 采样 样本 数 由 
459 增加 至 999 时 ， 频 谱 分 辩 率 由 0. 4357MHz 提升 至 0. 2002MHz， 基 于 MWC 采样 系统 的 总 采样 样本 
数 仅 增 加 至 222， 且 重 构 误差 在 10 悦 量 级 、 重 构 时间 在 10 习 量 级 。 图 2 中 黑色 频谱 是 低频 射电 天 文 
言 号 的 频谱 图 ， 红 色 频 谱 是 k 分 别 为 17、27、37 时 重 构 信号 的 频谱 图 ， 通 过 图 2 可 以 看 出 ， 基 于 
MWC 采样 系统 重 构 出 的 频谱 随 着 的 K 增加 ， 频 谱 分 辨 率 得 到 了 提升 ， 频 谱 带 宽 逐 渐 收 敛 并 接近 低 
频 射电 天 文 信号 。 因 此 ， 采 用 MWC 采样 系统 ， 不 仅 能 实现 单个 通道 低速 率 采样 ， 而 且 能 通过 增加 
时 域 采 样 时 长 来 提升 信号 的 频 域 分 辩 率 。 

#2 不 同 K 值 的 仿真 实验 数据 


m 


Tab. 2 The simulation experimental data at different K values 


单个 通道 采样 样本 数 17 27 37 
MWC 采样 系统 采样 样本 数 102 162 222 
等 效 总 采样 样本 数 459 729 999 
频谱 分 辩 率 /MHz 0. 4357 0. 2743 0. 2002 
重 构 误差 0. 035 0. 063 0. 028 
OMP 算法 的 重 构 时 间 /s 0. 0039 0. 0049 0. 0058 


T 


T T 
低频 射电 天 文 信号 
一 - 米 一 -- k 为 17 时 


T T 
低频 射电 天 文 信号 
=- k27 


T 


T T 
低频 射电 天 文 信号 
一 - 米 一 … k 为 37 时 


图 2 不 同 K 值 的 频谱 


Fig.2 The spectrum at different K values 


4 结论 
x 本 文通 过 对 MWC 采样 系统 的 深入 分 析 ， 以 及 实验 仿真 ,验证 了 增加 MVC 采样 系统 的 时 域 采 样 时 
2 长 ， 以 实现 低速 率 采 样 条 件 下 ， 提 升 信号 频 域 分 辨 率 的 可 行 性 和 有 效 性 。 为 提升 高 频 窄带 信号 和 低 
^ 频 射 电 天 文 信号 频谱 分 辩 率 提供 一 种 新 的 方法 和 途径 ,在 实际 工程 应 用 中 ， 需 结合 信号 特征 ， 进 一 
= 步 考虑 输入 信号 的 频带 带宽 、 频 谱 波形 、 奈 硅 斯 特 带宽 等 ， 设 计 WMC 采样 系统 的 通道 数 、 伪 随机 系 
© 列 周 期 、ADC 采样 周期 以 及 压缩 感知 矩阵 等 重要 参数 对 系统 性 能 的 影响 。 下 一 步 工作 将 针对 射电 文 

这 号 的 频谱 特征 和 频谱 分 辩 率 的 要 求 ， 设 计 MWC 采样 系统 ， 硬 件 实现 该 采样 系统 并 进行 实测 ， 力 人 


推进 论文 研究 成 果 的 实际 应 用 。 
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Research on Spectrum Resolution Based on Modulated Wideband 


Converter 
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(1. School of Information Science and Engineering, Yunnan University, Kunming 650500, 
China; 2. Yunnan Astronomy Radio Astronomy and VLBI Research Group, Chinese Academy of 
Sciences, Kunming 650216, China) 

Abstract: At present, the research on compressed sensing has mainly focused on reducing 


the reconstruction error and the reconstruction time by improving the reconstruction 


algorithm. However, the reconstruction error and reconstruction time lie in the 
practically allowed range. But researchers less pay attention to improving techniques of 
spectral resolution in the frequency domain analysis. In this paper, we use the modulated 
wideband converter (MWC) to sample the sparse signals, and employ the orthogonal matching 
pursuit (OMP) algorithm to reconstruct the input signals. By extending the sampling time 
to realize low sampling rate in a channel, and improve the spectral resolution. The 
exper imental results of low-frequency radio astronomical signals show that show that the 
spectral resolution increases from 0. 4357MHz to 0. 2002MHz when the number of samples in 
time domain increases from 459 to 999 with the same sampling rate for all channels. 
Meanwhile, the total sample number of the system increases only from 102 to 222. Our work 
provides a good reference for processing the narrowband high-frequency signal, such as 
radio astronomical signals. 

Keywords: Modulated Wideband Converter, Spectral Resolution, Compressed Sensing, 


Radio Astronomical Signal 


